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Usual Woods-Saxon single particle levels with BCS pairing are not able to reproduce the exper- 
imental occupation probabilities of the proton and neutron levels IP3/2, lpi/2, Of 5/2, 0g9/2 in the 
double-beta decay system '^^Ge to ^®Se. Shifting down the 0(79/2 level by hand can explain the data 
but it is not satisfactory. Here it is shown that a selfconsistent Hartree-Fock+BCS approach with 
experimental deformations for ^^Ge and ^''Se may decisively improve the agreement with the recent 
data on occupation probabilities by Schiffer et al. and Kay et al. Best agreement with available 
data on ^®Ge and '^^Se, as well as on neighbor isotopes, is obtained when the spin-orbit strength for 
neutrons is allowed to be larger than that for protons. The two-neutrino double-beta decay matrix 
element is also shown to agree with data. 

PACS numbers: 23.40.Hc, 21.60.Jz, 27.50.-|-e 

Among double-beta decay with two neutrino processes, the best studied case is that of ^^Ge going to ''^Se. Within 
a sphericahy symmetric description of these nuclei, the relevant proton and neutron single particle levels involved in 
the double-beta decay process are in the valence shells IP3/2, lpi/21 O/5/2, ^99/2- Their occupation probabilities have 
been recently measured by Schiffer et al. T] (for neutrons) and by Kay et al. [2^ (for protons). These data are not 
reproduced by previous mean field calculations using Woods-Saxon potential and pairing within BCS approximation. 
By increasing the binding energy of the Ogg /2 level the agreement is better, since its occupation probability increases 
and eventually reaches the experimental one. This energy shifting has been performed by hand in previous works 
[3-5 to obtain the experimental occupation probabilities, and the resulting ground state structure has been used 
to compute neutrinoless double-beta matrix elements. Shell-model calculations [6] have also been carried out with 
occupation probabilities in good agreement with experimental ones, but the single particle energies used are not 
available in the literature. 

In this work we use a deformed Skyrme Hartree-Fock (HP) mean field with pairing correlations within BCS approx- 
imation. At each HP iteration the BCS equations are solved to get new single particle wave functions $i, energies 
Cj and occupation probabilities until convergence is achieved. After convergence, the proton-neutron quasiparticle 
random phase approximation (pnQRPA) equations are solved on each deformed ground state basis for ^^Ge and ''^Se 
to get their Gamow- Teller (GT) strength distributions and to compute the two-neutrino double-beta decay matrix 
element. The formalism is similar to that used in previous works [7l-(9] with an important difference. In this work we 
allow the Skyrme force parameters to vary and search for a parametrization that provides better agreement with the 
experimental occupations of the above mentioned valence shells in ^^Ge and ^^Se, as well as in their neighbors. We 
focus on the effect of the spin-orbit strength, Wq, that we find to be of primary importance for the desired description 
of proton and neutron valence shells occupations. 

Experimental data and HP-fBCS calculations suggets prolate deformations for the ^^Ge and ^^Se ground states. 
The most precise data on the quadrupole moments of these nuclei come from Coulomb excitation reorientation 
measurements [10| which, together with experimental charge radii |11] . give a deformation parameter j3 =0.10 for 
^^Ge and /3 =0.16 for ^^Se (with 25-30% uncertainty). These experimental deformations are taken into account when 
obtaining our theoretical occupation probabilities in the active shell, and naturally leads us to check if the quadrupole 
deformation is responsible (fully or in part) for the disagreement between previous theoretical occupations and the 
experimental ones. To analyze this we show in Pig.[l]the experimental occupations in the active shell together with our 
theoretical results for different values of the quadrupole deformation, including the experimental and selfconsistent one, 
constraining the quadrupole moment in the HP-t-BCS procedure with Sk3 Skyrme force |12) . Since the single-particle 
orbital and total angular momenta are not good quantum numbers in the deformed cases, a projection to a spherical 
basis has been performed to obtain the occupation probabilities of the spherical valence shells. The uncertainties of 
the experimental data, not shown in the plot, are of a few percent (less than 10%), except for the highly unoccupied 
I — 4 proton state whose uncertainty is 30% in Se and over 100% in Ge. As the quadrupole deformation increases, 
the trend that can be observed in the figure is a decrease in the occupation probabilities of the states / = 1 (IP3/2 
and lpi/2 together) and I = 3 (O/5/2) and an increase in the occupation of the state I — 4 (0(79/2). Although this 
effect goes in the right direction, there are some features of the experimental occupations that are not reproduced by 
the deformed calculations as long as the deformation is constrained to reasonable values agreeing with experiment. 
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The most noticeable failure of the theory is a lack of occupancy in the 059/2 neutron level. We have checked that 
this failure is shared by a large variety of the Skyrme forces where the Skyrme spin-orbit strength Wq is equal for 
protons and neutrons. In any case, the selfconsistent deformed single nucleon potential gives better agreement with 
experiment for the occupation probabilities than a Woods-Saxon potential 'W, . 

To improve agreement with experimental occupations, we found it crucial to vary the spin-orbit terms of the Skyrme 
energy density functional. More specifically, to increase the value of the Wq parameter in the HF+BCS potential for 
neutrons. 

It has been already pointed out in the past that the Skyrme interactions do not reproduce properly the spin- 
orbit splittings of some nuclei, especially concerning neutron states, and that a more flexible two-body spin-orbit 
interaction is required |13j . The effect of modifying the Sk3 spin-orbit parameter Wq for neutrons is shown in Fig. [2j 
which contains results from a HF(Sk3'')-|-BCS calculation for spherical shape. The generalized Skyrme interaction 
SkB" maintains the original SkS parametrization except for the value of the spin-orbit strength Wq for neutrons that 
has been replaced by a new strength Wi, affecting only the terms of the potentials containing neutron densities. The 
value of the neutron strength has been modified from an initial value W,y = Wq =120 MeV fm^ all the way up to the 
value Wi, =240 MeV fm^. The corresponding strength for protons has been kept fixed, Wyr = Wq =120 MeV fm^ 
(the units of this parameter will be omitted in the following discussions). It is worth pointing out that because of 
self consistency the change of Wq for neutrons affects not only the neutron levels but also de proton levels. 

As shown in the plot, a good agreement on neutron occupation probabilities between theory and experiment is 
found for Wi, =200; in this situation the Ogg/2 state is more bound than the lpi/2 state. The result for protons yields 
also a good agreement with experiment; in this case, the proton state lpi/2 is still more bound than the state Ogg/2- 

In Fig. [3] one can see the proton and neutron single particle energies in the active shell corresponding to the 
occupation probabilities of Fig. [2] together with the Fermi levels. As the spin-orbit strength for neutrons increases, 
the single proton and neutron levels lpi/2 and O/5/2 get clearly less bound, whereas the level 0g()/2 gets more bound. 
The level 039/2 exchanges its energy position with the level lpi/2 in the vicinity of W„ =240 in the case of protons 
and for 160< <200 in the case of neutrons. In this last case, the level 0^9/2 also exchanges its energy position 
with the level O/5/2 but in the range 200< Wi, <240. Finally, the level IP3/2, which is in general the most bound of 
them, remains with about the same binding energy. In what follows we refer to the interaction with Wti — 120 and 
= 200 as Sk3^. 

The considerable improvement obtained by setting different values of the spin-orbit interaction for each type of 
nucleon supports a modification of the spin-orbit term in the general Skyrme interaction 14 , that could have a 
structure similar to that of the other Skyrme terms, namely containing a factor of the form (1 + eP^.), where Pr is an 
isospin-exchange operator. More involved isospin-dependent formulations are also possible and worth exploring. 

We show in Fig. |4]the occupation probabilities in the active shells obtained with spherical IIF(Sk3)-('BCS and with 
the improved IIF(Sk3'^)-|-BCS (using VF^ =120 and W^, =200) and the experimental deformation. The results are also 
compared to the experimental data. All of them are normalized so that the total number of nucleons in the active shells 
is the expected one: 4 protons and 16 neutrons in ^^Ge and 6 protons and 14 neutrons in ^^Sc. In our calculations 
the actual numbers of protons and neutrons in these valence shells are lower, especially in the deformed cases, where 
the occupation probabilities are spread along many spherical single particle states, not only in those active shells 
of the spherical shell model. When obtaining the occupations from experiment, a similar (although not identical) 
renormalization procedure was also performed (see [l][2] for details). As seen in the figure, the improvement reached 
by combining the effect of deformation and modified spin-orbit strength with respect to the spherical HF(Sk3'^)-|-BCS 
is apparent. It is possible however that other ingredients are still missing which could enhance the trend initiated 
by deformation and spin-orbit strength towards a reliable description of the single particle spectrum in this nuclear 
region. 

This new spin-orbit strength can be applied to the neighbor nuclei ''^Ge and ''^Se, whose active-shell occupation 
probabilities have been also measured [II [2] . In Fig. [s] we compare the measured probabilities with our theoretical 
results for the IIF(Sk3)-|-BCS calculation with spherical shape and for the new HF(Sk3*)-|-BCS calculation with 
experimental deformations. The /3-values extracted from data in refs. TO', TT are /3 = 0.13 for ^"'Ge and f3 = 0.12 for 
^®Se. The inclusion of experimental deformation and a larger spin-orbit strength for neutrons change also in this case 
the theoretical results towards a better agreement with the experimental data. 

The new single particle energies and occupation probabilities have also an effect on the matrix element of the 
two-neutrino double-beta decay ""^Ge— >^^Se. In Fig. [g] we show the theoretical double-beta decay matrix element 
running sum as a function of the excitation energy of the intermediate nucleus ^^As. The results shown correspond to 
deformed pnQRPA calculations (see ref. [S] and references therein for details of calculations) using the single particle 
basis for ^^Ge and for ^^Se obtained with the new {SkS", > VF^) and the standard (Sk3, = W„ = Wq) 
parametrizations, with and without deformation. Proton-neutron particle-hole ph and particle-particle pp residual 
interactions within QRPA are included to obtain the contributions to the double-beta matrix element, following the 
procedure explained in previous works [S] . The coupling constants of the ph and pp channels take the values Xph — 
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0.25 MeV and Kpp = 0.08 MeV respectively. In the case of Xph, it is an average value between the one used in previous 
works [H HSl [16] and the selfconsistent value [17]. In the case of Kpp, the current value is in agreement with the 
one used in some previous works [HI dS] , but it is twice as large as the one in [3] . We note that the larger value of 
the pp residual interaction coupling constant is responsible for the particular shape of the plots in Fig. [6] where the 
matrix element increases in the region from to 4 MeV of excitation energy, then goes up and down with increasing 
energy and finally decreases until reaching its final value. This behavior does not show up for the smaller values of 
the pp coupling constant in [S]. This finding also agrees with the conclusions reached in ref. |18j . namely that the 
characteristic increase and decrease with energy of the matrix element found in many shell model calculations is also 
found in our deformed pnQRPA calculations when we use larger values of Upp than those used in our previous works 

The experimental data shown in Fig. |6] arc obtained in a phenomenological way using the experimental Gamow- 
Teller (GT) strengths of the transitions ^Ge->'^^As (GT") [TS] and '^'^Se^'^^As (GT+) [ID], as explained in previous 
experimental or theoretical works [3 [2D]. An experimental range for the final value of the matrix element is also 
shown in the plot, obtained from the experimental double-beta decay half-life [21] with two possible values of the 
axial-to-vector ratio gA = 1.25 and = 1.00 (quenched value). 

It is clearly seen in the figure that the calculations with the increased spin-orbit strength for neutrons reproduce 
much better the experimental values. In the low energy region the agreement is improved with deformation. There are 
no large differences between the calculations with spherical ground states and the one with experimental deformations. 
For SkS"^ both of them lie within the experimental range for the total value of the matrix element. For Sk3 a large 
disagreement is found with the Kpp value used here. 

We can trace back the origin of the low-energy behavior of the double-beta matrix element, studying the accumulated 
strength distributions of the GT~ transition ^^Ge— )-''^As and of the GT+ transition ^^Se— )-^^As. In the case of the 
GT~ strength distribution, our result with Wi, > W-^ is roughly twice as large as the result with the usual strength 
(Wy = Wtt = Wo) until 4 MeV, and the former reproduces better the experimental data in this range. In the case 
of the GT+ strength distribution, the results with increased spin-orbit strength are again larger, and there is a big 
difference between the spherical or the deformed single particle basis for ^^Se. The deformed one gives the best 
agreement with the results using experimental data. Beyond 2 MeV the experimental value increases, whereas the 
theoretical results wait until 5 MeV to do so (not shown in the plot). 

From this analysis we conclude that the better agreement between theoretical and experimental occupation prob- 
abilities in the valence states with deformation and enhanced neutron spin-orbit gives rise to GT transitions in good 
agreement with the experimentally known strengths at low energy that involve the valence shells. This in turn produce 
a better agreement for the double-beta matrix element. 

Apart from the single- and double-beta decay properties of the nuclei under study, the modification of the neutron 
spin-orbit strength affects some bulk properties of the nuclei, in particular the binding energy and the charge radius. 
Since the parameters of the Skyrme interactions are fitted to reproduce these bulk properties in a given set of nuclei, 
the modification of one of them without an overall refitting of the rest may result in unrealistic values of those bulk 
properties. Notwithstanding the above, when comparing with the experimental binding energies of 661.6 MeV for 
^^Ge and 662.07 MeV for ^^Se we get a 0.5% underestimation in the calculation with the usual Sk3 force and a 2% 
overestimation in the calculation with the Sk3'' force. Concerning the charge radii, the IIF(Sk3)-|-BCS results are 
4.12 fm for ^^Ge and 4.17 fm for ''^Se, whereas with increased spin-orbit strength one gets 4.09 fm for ^^Ge and 4.14 
fm for ''^Se, which actually agree better with the experimental values of 4.08 fm for ^^Ge and of 4.14 fm for ^^Se [11) . 

To conclude, we have shown that agreement with experimental occupation numbers in the valence shells of ''^Ge and 
^^Se can be obtained within the selfconsistent Skyrme-t-HF-l-BCS deformed mean field approximation. Best agreement 
with experiment is obtained when a spin-orbit interaction more flexible than the standard one in Skyrme forces is 
considered, allowing in particular for a larger neutron spin-orbit force. We point out that the interplay between 
spin-orbit strengths and deformation is an essential feature to give the experimental occupations. These occupation 
numbers are in turn of great importance for beta strength distributions and two-neutrino double-beta decay matrix 
element calculations, especially in the low energy region. Changes at the level of the two-body Skyrme interaction 
and energy density functional are suggested to further improve the theoretical description of nuclei and will be further 
investigated. 
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FIG. 1: (Color online) Active-shell experimental [TJ |2] and theoretical occupations for different quadrupole deformations 
(indicated by the deformation parameter /3), from a HF(Sk3)-|-BCS calculation in ^"Ge and ^^Se. Occupations are shown for 
the proton and neutron states in the active shells, which include the single particle states lpi/2 and IP3/2 (gathered as I — 1), 
0/5/2 (/ = 3) and O59/2 (^ = 4). 
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FIG. 2: (Color online) Active-shell experimental [T][2] and theoretical occupations for different values of the neutron spin-orbit 
strength W^, (with VF^ =120), from a HF(Sk3'')-|-BCS calculation in ^^Ge and '^^Se with spherical ground states. 
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FIG. 3: (Color online) Theoretical Fermi levels and active-shell single particle energies for different values of the neutron spin- 
orbit strength (with =120), from a HF(Sk3'')-|-BCS calculation in ''''Ge and ^^Se with spherical ground states. The 
occupations corresponding to these levels appeared in Fig. |2] 
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FIG. 4: (Color online) Active-shell experimental [HIS] and theoretical occupations in ^"Ge and ^®Se. The first theoretical result 
is the original HF(Sk3)-|-BCS calculation for spherical ground states and the second one is the HF(Sk3'')-|-BCS calculation 
with Wi, —200 and experimental ground-state deformations [T0|. Both experiment and theory have been normalized so that 
the active shell contains all the nucleons exceeding the magic number 28. 
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FIG. 5: (Color online) Same as in Fig.|4]but for ^''Ge and ^®Se. 
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FIG. 6; (Color online) Running sum of the matrix element for the two-neutrino double-beta decay ^^Ge— >^®Se, as a function of 
the intermediate excitation energy in '^^As. The single particle basis for ^®Ge and ^''Se have been obtained from HF(Sk3'')-|-BCS 
{Wi, > Wtt) and from HF(Sk3)-|-BCS {Wv ~ W-n = Wo) calculations for spherical (sph) and experimental deformed (exp def) 
ground state shapes. An experimental range of the total matrix element is shown, coming from the experimental half-life of 
the process. The experimental low-energy steps of the matrix element running sum come from experimental GT^ and GT~ 
strengths jSD]. 



